
APPROXIMATION FOR 

RADIANT ENERGY TRANSPORT 

IN NONGRAY, NONSCATTERING GASES 



Lewis Research Center 
Cleveland, Ohio 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


WASHINGTON, D. C. 


JUNE 1967 


TECH LIBRARY KAFB, NM 



TECH LIBRARY KAFB, NM 



0130832 


APPROXIMATION FOR RADIANT ENERGY TRANSPORT IN 

NONGRAY, NONSCATTERING GASES 

By R. W. Patch 

Lewis Research Center 
Cleveland, Ohio 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


For sale by the Clearinghouse for Federal Scientific and Technical Information 
Springfield, Virginia 22151 — CFSTI price $3.00 



APPROXIMATION FOR RADIANT ENERGY TRANSPORT IN 
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SUMMARY 

Radiant heat transfer is the most important mode of energy transfer in many high- 
temperature gases and plasmas; thus, it is necessary to be able to estimate its magni- 
tude with confidence. In cases where the temperature varies appreciably over a photon 
mean free path, a transport analysis is generally required. When a plasma has a com- 
plex spectrum, such an analysis can be unduly time consuming even with a high-speed 
digital computer unless approximations are made. 

An approximation for integrated intensity in nongray, nonscattering gases is pre- 
sented that uses a new mean absorption coefficient in the gray-gas transport equation. 
This absorption coefficient is called an "effective absorption coefficient. " The approxi- 
mation was tested by calculating the integrated intensities by exact and by approximate 
methods for 39 cases. The integrated intensities calculated by the effective absorption 
coefficient approximation had errors of only -25 to +28 percent, which demonstrated the 
superiority of this approximation over a similar approximation using the Planck mean 
absorption coefficient. Despite certain restrictions on the form of the spectral absorp- 
tion coefficient, the effective absorption coefficient approximation is expected to.be use- 
ful for most problems of nonscattering plasmas in local thermodynamic equilibrium 
(LTE) even where large temperature and pressure differences and inhomogeneities are 
present. 


INTRODUCTION 

In very -high-temperature propulsion devices such as gaseous nuclear rockets, the 
dominant mechanism of heat transfer is radiant energy exchange between volumes of 
plasma and between plasma and the wall. There are two general methods of analyzing 
such heat transfer: transport analysis and the diffusion approximation. An introduction 
to these two methods is contained in reference 1. 



The term ’’transport analysis” generally refers to a method of calculating radiant 
flux density by considering conditions at other points in space besides the point at which 
the radiant flux density is desired. In this report, transport analysis refers to single or 
multiple use of the transport equation (eq. (131) of ref. 1 or some modification of it) to 
calculate integrated intensity in some particular direction in space from the point where 
the radiant flux density is desired. This calculation is repeated for other directions, 
and an integration over solid angle is performed to obtain the radiant flux density. If no 
approximations are made, transport analysis gives the correct value of radiant flux den- 
sity. 

The usual diffusion approximation is really an approximate transport analysis in 
which the radiant flux density is inversely proportional to the Rosseland mean absorption 
coefficient and directly proportional to the temperature gradient at the point where the 
radiant flux density is desired (ref. 1). This approximation cannot be equaled for sim- 
plicity, but it is not applicable when temperature varies appreciably in a photon mean 
free path as it may far from the center of the core of a gaseous nuclear rocket. 

The gray-gas transport equation is merely the monochromatic transport equation 
(eq. (131) of ref. 1) integrated over all wave numbers with the assumption that the spec- 
tral (or monochromatic) absorption coefficient is independent of wave number. The 
gray-gas transport equation is, in general, rigorously applicable only to a gray gas or 
to a gas that is optically thin at all important wave numbers (important wave numbers 
are those within about one order of magnitude of the wave number of the maximum spe- 
cific intensity of the Planck function) . 

For other cases, the monochromatic transport equation is customarily applied at a 
sufficiently large number of wave numbers to enable calculation of integrated intensity by 
numerical quadrature. The radiant flux density is then obtained by integrating the inte- 
grated intensity over all solid angles. If the absorption coefficient has an intricate de- 
pendence on wave number and the problem does not involve slab geometry, calculation of 
radiant flux density by this method may be too time consuming even with a high-speed 
digital computer. For a uranium plasma, this difficulty can be appreciated as more than 
9000 spectral lines have been observed in a uranium arc. For this and other reasons, 
radiant heat transfer far from the center of the core of a gaseous nuclear rocket has not 
been analyzed with any confidence. 

Attempts to find satisfactory approximations for calculating radiant energy transport 
in nongray gases of any optical thickness have been made by a number of investigators. 
Stewart (ref. 2) has obtained a solution for slab -geometry steady-state problems where 
the spectral absorption coefficient can be factored into a part that is dependent on tem- 
perature and density and a part that is dependent on wave number. Grant (ref. 3) has 
obtained an approximation by dividing the spectrum into wave -number (or frequency) 
groups and by dividing the spectral absorption coefficient into absorption coefficient 
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groups. Unfortunately, the latter approximation is only valid when the spectral absorp- 
tion coefficient is independent of position. Sampson (ref. 4) has obtained an approxima- 
tion for slab-geometry problems that gives radiant flux densities and their divergences 
that are * 'nearly always within a factor of 2 of the correct values. " This approximation 
gives exact values of these quantities for three classes of problems: (1) an optically 
thick gas where the diffusion approximation is applicable, (2) a gas that is optically thin 
at all important wave numbers, and (3) a gray gas. All three methods lack either gener- 
ality or accuracy. Viskanta (ref. 5) and Campbell and Nelson (ref. 6) used flux -weighted 
mean absorption coefficients that require calculation of spectral fluxes in order to evalu- 
ate them. 

In this report, a new approximation is presented in which the integrated intensity is 
obtained from the gray-gas transport equation by using therein an effective absorption 
coefficient averaged over the spectrum once and for all. The use of the effective ab- 
sorption coefficient in the gray-gas transport equation is expected to give a satisfactory 
approximation to integrated intensity for most cases involving a nonscattering gas in 
local thermodynamic equilibrium and possessing spectral lines, continua, or both. The 
effective absorption coefficient gives exact values of integrated intensity for five classes 
of problems. The effective absorption coefficient approximation was successfully ap- 
plied to calculate integrated intensity for 39 cases, 33 of which did not belong to the five 
classes. 


SYMBOLS 

A area, cm 

Aa perturbation in spectral linear absorption coefficient, cm”* 

a g effective linear absorption coefficient, cm”* 

a_ gray linear absorption coefficient, cm"* 

S _i 

a pl Planck mean linear absorption coefficient, cm 
a Ro Rosseland mean linear absorption coefficient, cm”* 

ay spectral (or monochromatic) linear absorption coefficient at photon energy V, 
cm”* 

a~ spectral (or monochromatic) linear absorption coefficient at wave number v, 

v -1 
cm 

B integrated Planck function (called integrated blackbody intensity in ref. 1), 
W/(cm ^ ) ( ste radian) 

B~ Planck (or blackbody) function, W/(cm)(steradian) 
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col mi 


b collision half -width of Lorentz line (line with a dispersion contour), cm' 1 
b average collision half -width of Lorentz lines, cm”* 
c velocity of light in a vacuum, cm/sec 

o 

E~ power radiated by spontaneous emission, W/(cm )(steradian) 

— ► o 

F radiant flux density, W/cm 

F n component of F in the direction of n, W/cm 

f cumulative distribution 

f~* inverse function of f 

h Planck constant, (J)(sec) 

o 

I integrated intensity, W/(cm )(steradian) 

I~ specific intensity, W/(cm)(steradian) 

J line strength parameter 

k Boltzmann constant, J/°K 

L distance from point P to wall (fig. 1), cm 

M parameter in series expansion of Q 

n unit vector normal to dA (fig. 1) 

2 

P probability parameter, cm 

p pressure, atm 

Q fraction of Planck function between V and infinite v 

q integral in denominator of effective absorption coefficient for narrow wave-number 

group, cm”* 

R random number 

S distance along ray (fig. 1), cm 

_2 

integrated linear absorption coefficient of spectral line, cm 

-2 

average integrated linear absorption coefficient of spectral lines, cm 
T temperature, °K 

t distance- or temperature -dependent factor of spectral linear absorption coefficient 
V dimensionless photon energy 

W square of V 
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w wave-number-dependent factor of spectral linear absorption coefficient, cm' 

Aw perturbation in w, cm - * 

w g constant term in wave-number-dependent factor of spectral linear absorption 
coefficient, cm - * 

X mole fraction 

x component of S in direction n (fig. 1), cm 

y dummy variable 

z integral in numerator of effective absorption coefficient for narrow wave- 

number group, cm - ^ 

6 angle between ray and direction n (fig. 1), rad 

V wave number, cm - * 

Av difference of two wave numbers, cm - * 

a Stefan-Boltzmann constant, W/(cm^)(°K)^ 

r g effective optical depth 

*7* gray optical depth 

r~ spectral (or monochromatic) optical depth 

ip local conditions in gas (temperature and pressure or, if more than one gas, 

partial pressures of the different gases) 

solid angle (fig. 1), steradians 

Subscripts: 

a gas a 

approx approximate 

b gas b 

d gas d 

e gas e 

i isothermal, homogeneous, and constant pressure 

j alternate subscript for origin group or gas 

k alternate subscript for origin group or gas 

l alternate subscript for origin group or gas 

line spectral line 
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m wave -number group, minimum wave number in wave -number group m, or mini- 

mum photon energy in wave -number group m 

m+1 maximum wave number in wave -number group m or maximum photon energy in 
wave -number group m 

max maximum value 

min minimum value 

n component in direction n (fig. 1) 

O s * 0 

o center of spectral line 

r origin group, species, or gas 

u number of origin groups or gases 
v alternate subscript for origin group or gas 
w wall 

a number of wave -number groups 


ANALYSIS 

The derivation of the effective absorption coefficient for a single gas or plasma, to- 
gether with refinements for large temperature ratios and for inhomogeneous mixtures, 
is contained in this section. The effective absorption coefficient is derived initially for 
use in the gray-gas transport equation in order to express the integrated intensity for a 
homogeneous, isothermal, isobaric, nongray gas. A hypothesis is then made that allows 
the use of this equation as an approximation to the integrated intensity for a gas with any 
kind of gradients. 

This hypothesis has two principal defects, one of which causes appreciable errors 
for large temperature ratios and the other causes errors for inhomogeneous mixtures. 
The first defect is largely remedied by dividing the spectrum into wave -number groups, 
with a modified effective absorption coefficient and modified gray-gas transport equation 
for each group. The second defect is largely remedied by dividing the spectral absorp- 
tion coefficient into origin groups according to the origin of the spectra and by assuming 
no correlation between the spectra of different origin groups. This division requires the 
assignment of an effective absorption coefficient to each origin group. 

The foregoing approximations and hypothesis are extended to take into account a 
black wall at finite temperature. 
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Assumptions and Restrictions 

The assumptions utilized in the following derivations are listed herein for conven- 
ience. 

(1) The term "gas" includes ’’plasma." 

(2) No scattering of radiation occurs. 

(3) The gas is in local thermodynamic equilibrium, that is, the composition, degree 
of molecular excitation, and velocity distribution at any point correspond to thermo- 
dynamic equilibrium at some temperature, but the specific intensities do not have to 
correspond to the Planck function at that temperature. 

(4) The gas is nonrelativistic. 

(5) Changes in the state of the gas are slow enough that the time of flight of photons 
may be neglected. 

(6) Doppler shift in the spectrum may be neglected. When narrow spectral lines are 
present, this is a more stringent requirement than (4) or (5). 

(7) Any walls are black, that is, they do not reflect or transmit light. 

(8) The ordinary index of refraction is 1. 

(9) There is no interaction between photons. 

(10) Only one photon interacts with a particle at any given time. 

(11) The molecules of the gas have random orientation in space. 

The analysis is neither limited to gases with only spectral lines or only continua, nor is 
it restricted as to optical thickness. 


Basic Relations 

Some fundamental relations are needed in subsequent sections and are given here. 
Let a~ be the spectral linear absorption coefficient of a gas, including the induced 
emission factor (ref. 7). If S is the distance along a ray from the point P (fig. 1), 
then the spectral optical depth of S relative to P is 

rs 

T v S J 0 a ? (S ’ )dS ' (1) 

The Planck function is given by 

B~ = — 2hc2y3 (2) 

* e hc?/kT _ j 
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Here h is the Planck constant, c is the velocity of light, V is the wave number, k is 
the Boltzmann constant, and T is the temperature. 

The value of the specific intensity I~ at point P in the direction from the wall 
at L toward P is obtained by integration of equation (116) of Chandrasekhar (ref, 1) 
with appropriate boundary conditions. 

/-L -t~(S) -t~(L) 

% = J Q B~(S)e a~(S)dS + B~ w e (3) 

This is the monochromatic transport equation. The quantity B~ is the Planck function 

for the gas temperature, whereas B~ is the Planck function for the wall temperature. 

w 

The integrated intensity I at point P is found by integrating equation (3) over all 
wave numbers 

/ °° /'L -T~(S) /•«> -T~(L) 

J f B~(S)e " a~(S)dSd?+^ B~ w e ^ dV (4) 

According to equation (4), the contribution to I at point P due to the gas along an ele- 
ment of ray dS is 


dl = 




- T f?(S) 


a~(S)d?dS 


(5) 


For a gray gas the spectral linear absorption coefficient including the induced emis- 
sion factor (ref. 7) is independent of wave number and is designated a g . The corre- 
sponding gray optical depth at S is 



( 6 ) 


Hence, for a gray gas, equation (4) becomes 



B(S)e 


■t (S) 

g a g (S)dS 


B e 
w 


-UL) 


(7) 


where B is the integrated Planck function for the gas temperature, and B w is the inte- 
grated Planck function for the wall temperature. According to equation (7), the contribu- 
tion to I at point P due to the gas along an element of ray dS is 


-t (S) 

dl = B(S)e g a (S)dS = B(S)e 

o 


-^V' )dS ' 


a g(S)dS 


( 8 ) 
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Effective Absorption Coefficient for Single Gas 

In this section, the effective absorption coefficient is derived for an isothermal, iso- 

baric, single gas or homogeneous mixture, and its mathematical behavior is discussed. 

A hypothesis is made that permits its application to gases with any kinds of gradients. 

Homogeneous, isothe rmal, isobaric, nongray gases . - In order to use equations of 

the form of equations (7) and (8) for a homogeneous, isothermal, isobaric, nongray gas, 

an effective absorption coefficient a must be found to use in place of a . Specifically, 

e g 

-/■ Sl a e (s;)ds; 

dl = Be 0 a e (S.)dS. (9) 

where the subscript i has been added to S to indicate that the gas is homogeneous, iso- 
thermal, and isobaric. The effective absorption coefficient a g must take whatever form 
is necessary so that equation (9) is valid for any homogeneous, isothermal, isobaric gas. 
For a given wave number, a~ and B~ are independent of S^; thus equating equations 
(5) and (9) results in 



V V 


-a~S. 
v 1 


dv = Be 


-jfVlK 


,< s i> 


Integrating equation (10) with respect to from 0 to sV gives 



B~e 


. a <vS" 

V l 


dV + B = -Be 


/'S! 

-J 0 ‘>,tM 


+ B 


(10) 


( 11 ) 


Cancelling the two B's, dropping the double primes in equation (11), and dividing equa- 
tion (10) by equation (11) yield 


a e (S t ) = 


X 

X 


-a~S. 

B^a^e v 1 dV 

oo -a^S. 

B^e v 1 d?^ 


( 12 ) 


For a homogeneous, isothermal, isobaric, gray or nongray gas and a cold wall at L, 
equation (9) can be integrated to give 
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( 13 ) 



where the effective optical depth r g is defined by 

T e< S i> * f‘ a e< s l) <ls l < l4 > 

Substituting equation (12) into equation (14) and noting that the numerator in equation (12) 
is just the negative of the derivative of the denominator with respect to Sj give 


r e (S.) = in 


B 



B~e 


~ a u S i 


d? 


(15) 


The effective optical depth T g is always measured from the point P where the inte- 
grated intensity and radiant flux density are desired. 

Behavior of a g for a homogeneous, isothermal, isobaric gas. - The variation of 

a e as ®i varies I s relatively simple. As approaches 0, it is clear from equa- 
tion (12) that a g approaches the Planck mean absorption coefficient ap^ given by 


api ‘lX Bv&v ** (l6) 

As Sj approaches infinity, a g approaches the smallest value of a~ in the spectrum. 
For intermediate values of the function a~ which results in a maximum value of a g 
with apj held constant can be shown to be the constant value a pl . Hence a g s apj 
always. 

The variation of a g with has a physical significance for a nongray gas. When 
Sj is very small, virtually all the radiation emitted by gas along dS i will reach point P 
(fig. 1). The appropriate mean absorption coefficient for use in the gray-gas transport 
equation for such optically thin cases is known to be apj (ref. 4) because there is no 
measurable self -absorption, that is, the gas does not absorb its own radiation to any 
extent. 

If the gas is gray and if r g is large, then radiation emitted by gas along dSj will 
be strongly self-absorbed before it reaches point P; thus, only a small fraction of the 
radiation emitted along dSj gets to P. All photons emitted along dSp regardless of 
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wave number, will have the same probability of getting to P, namely, exp(-Tg). Equa- 
tion (7), which contains exp(-T ) in its first term, then gives the correct integrated in- 
tensity at P. 

If, on the other hand, the gas is nongray, the dependence of the emission of the gas on 
a^ must be considered. If equation (119) in reference 1 is used, the Kirchhoff law may 
be expressed as E~= a~B~ where E~ is the power per unit volume per unit solid angle 
per unit wave number radiated by spontaneous emission. In addition to the dependence of 
E^. on B^ E^. is large at those wave numbers where a~ is large. The absorption of 
radiation from a nongray gas by the same gas can therefore be termed "preferential 
self -absorption" because for a given path length the gas will absorb a larger fraction of 
its own radiation than of radiation from a blackbody at its temperature. 

For an optically thick, nongray gas, the contribution to the integrated intensity at P 
by gas along dS^ far from P is smaller than would be calculated if ap^ were used in 
equation (7) or (13). For a homogeneous, isothermal, isobaric, nongray gas, this effect 
due to preferential self -absorption is allowed for properly by using a g in equation (13) 
instead of apj because a g is weighted properly to account for preferential self- 
absorption. 

Nongray gases with gradients. - In cases where a~ can be approximately or exactly 
factored into a product w(V)t(S), the concept of preferential self -absorption is valid. 

The degree of preferential self -absorption (or, equivalently, how much less a g is than 
apj) would be expected to depend on the effective optical depth r g rather than S. 

The hypothesis is made here that r g can be used as an approximate index of the de- 
gree of preferential self -absorption in any gas of a single species whether or not a~ 
can be exactly factored into a product w(ty)t(S) and whether or not temperature and pres- 
sure are constant. The basis for this hypothesis is that the spectra of an actual gas at 
moderately different temperatures and pressures are not markedly different from one 
another except in magnitude; thus, in a practical problem a^ does approximately factor 
into a product w(£")t(S). This hypothesis is used below in finding an expression for a g , 
whose value relative to apj is merely a measure of the degree of preferential self- 
absorption in the gas. 

The value of a. may be computed if the value of a at S in a gas with gradients 
is assumed to be the same as the value of a g for an isothermal, isobaric gas having the 
same effective optical depth (same degree of preferential self -absorption) and same 
pressure and temperature as at S in the gas with gradients. This assumption may be 
called the equal-effective-optical-depth hypothesis. Before expressing this hypothesis 
mathematically, the quantity Sj is redefined as the distance along a ray in an isothermal, 
isobaric gas having temperature and pressure the same as at S and having such a value 
that the effective optical depths r in the gas with gradients and in the isothermal, iso- 
baric gas are the same. Then a g for a gas with gradients is given by the expression 
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for the associated isothermal, isobaric gas. 


a. 




/ oo 

/ oo 






a^)S. 

e v 1 dT* 


( 17 ) 


where xj/ stands for the temperature and pressure at S in the gas with gradients and, 
of course, everywhere in the associated isothermal, isobaric gas. 

The effective optical depth for a gas with gradients is defined by 

T e (S) = jT S a e [*(S’), S^S'jjdS’ (18) 

which is analogous to equation (14). The effective optical depth for an isothermal, iso- 
baric gas of depth with the same pressure and temperature as at S in a gas with 
gradients is, according to equation (15), 


T 


e 


= In 


B(S) 



(19) 


The definition of given previous to equation (17) may then be expressed mathemati- 
cally by equating equations (18) and (19). 



a e [^(S'), S.(S , )]dS’ h 


t = In 
e 


/' 

Jo 


B(S) 

-a~(S)S.(S) 
B~(S)e v 1 dV 


( 20 ) 


hi general is not equal to S. Solving equation (20) for S^(S) and then equation (17) 
for a (S) is a time consuming procedure by any method. This procedure may be avoided 

V 

by noting that, for a given value of S, S.^ and T g are uniquely related as shown by the 
right-hand side of equation (20). Hence a g may be regarded as a function of xp(S) and 
r g ; thus, differentiating the left-hand side of equation (20) gives 

= a e[*<*>> T e] < 21 > 
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and the quantity Sj no longer enters explicitly into the calculation of T g . Equation (21) 
may be integrated readily by numerical methods to get r e (S) if a e [i//(S), tJ is available 
in a form such as a table. This integration is started by using the boundary condition 
T g = 0 at S = 0. 

For a given gas, values of a e |i//(S), tJ may be tabulated once and for all in a three- 
dimensional table with arguments of pressure p, T, and T g because a 0 depends only 
on p, T, and r g . Such a table may be constructed by selecting p, T, and a number of 
values of S.. The quantities a. and t are calculated from equations (17) and (19), 
respectively, for each value of S., and the results are interpolated to find a at the de- 
sired values of r g for the table. This procedure is then repeated for other values of 
p and T. 

The previous considerations make it possible to generalize equation (13) to give the 
effective transport equation for a gas with or without temperature and pressure gradients 
and with a cold wall at L, namely, 



B(S)e 


• r;]dS’ 




,<U 


Be 


dr. 


( 22 ) 


When I must be evaluated a sizable number of times for the same gas and when that 
gas possesses a complex spectrum, it is faster to tabulate a g ^(S), rj once and for all 
and then to solve equations (21) and (22) once for each I than to evaluate the double in- 
tegral in the exact equation (4) once for each I. This advantage becomes clearer if it is 
considered that for a plasma with, say, 9000 spectral lines the integral with respect to S 
in equation (4) would probably have to be evaluated by numerical methods at 40 000 or 
more different wave numbers so that the integral with respect to wave number could be 
evaluated by numerical methods with reasonable accuracy. This procedure would have 
to be followed for each I if the exact equation (4) were used. 

Exact classes of problems . - Equation (22) has already been shown to be exact for 
(1) isothermal, isobaric, single gases and homogeneous mixtures and is shown in the 
next two paragraphs to be exact for (2) gray or (3) optically thin gases with or without 
temperature or pressure gradients and with a cold, black wall behind them. 

In the case of a gray gas, equation (17) gives a. = a_. Substituting this equality into 

“ o 

equation (22) and utilizing equation (6) gives equation (7) except for the last term, which 
is appreciable only if the wall is not cold. Since equation (7) is exact, it follows that 
equation (22) is exact for gray gases with a cold, black wall behind them. 

For a gas that is optically thin at all wave numbers, that is, a^S « 1» a e - a pi’ ® 
this substitution is made and the exponential is expanded in a power series, equation (22) 
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becomes 


= B(S) |l - f* ap^S’JdS* + . . .J j 


PI 


(S)dS 


Because a~S « 1 at all wave numbers, 



ap^S’JdS’ « 1 


so all but the first term of the power series may be neglected, giving 


I = 



B(S)a pl (S)dS 


(23) 


(24) 


(25) 


Some manipulation of exact equation (4) is necessary to show that it is identical to equa- 
tion (25) for an optically thin gas. If the wall is cold, the last term of equation (4) may 
be dropped. Setting the exponential equal to 1 in equation (4), reversing the order of in- 
tegration, and utilizing equation (16) gives equation (25). Hence, equation (22) is exact 
for a gas that is optically thin at all wave numbers and has a cold, black wall behind it. 


Refinement by Using Wave-Number Groups 

Error due to Planck function shift . - In problems involving large maximum -to- 
rn in imum temperature ratios, nongray gas, and spectral optical depths of order of mag- 
nitude 1 or more, an appreciable error is frequently present in equation (22) because of 
the shift of the curve of B~ against V to larger values of V as temperature increases. 
This error is not evident from a cursory inspection of equation (22) because the depen- 
dence on the Planck function B~ is hidden in the evaluation of a g (eq. (17)). The error 
can be demonstrated, but the derivation is long. The error is principally due to the inte- 
grand of the exponential of equation (22) being calculated with values of a g determined 
with Planck functions corresponding to temperatures at S* rather than at S. 

Introduction of wave-number groups . - The error may be minimized by writing sev- 
eral equations analogous to equation (22) but with each applicable to a particular wave- 
number range. Thus, every important wave number is included in one of the wave- 
number ranges, which will hereafter be called "wave-number groups. " The value of I 
may be obtained by simply adding the several equations. There is a precedent (ref. 3) 
for such a procedure. From the derivation that follows, it will become evident that when 
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I 


this procedure is carried to the extreme of having an infinite number of equations for an 
infinite number of very narrow wave -number groups, the value of I obtained is always 
exact. 

The group intensity I m of the m wave-number group is defined as all intensity 
with wave numbers between V m and and with the help of equation (3) ^ is given 

exactly by 


m 


s f I?d? = f f Tv 

Jo 


(L) 


B^( T ^) e 


- T ? 




dP + 


j£ b F, 


-r v (h) 


w' 


dV 


(26) 


where m under the integral signs means to integrate from V m to ^ 

If the same derivation is followed as for equation (17), except that 1^ is used 
wherever I was used, the group effective absorption coefficient is found to be 


a e,mH’ S l,m] = 


/ -a~(»^)S. „ 

v 1>m d? 

ill 

/ -a~(^)S. „ 
B~(^)e v 1>m d? 


(27) 


where an equal-effective-optical-depth hypothesis is required for each group. Thus, for 
the same value of S, values of are generally different for different groups. 

Several other group quantities result from such a derivation. The group Planck 
function is defined by 


m 


* f D v 

m 


dV 


(28) 


The group effective optical depth t q m is given by 


f Q S a e,mh s ’>’ s l,m< s ’>] ds ' * T e,m = 


B m< s > 


f -a~(S)S. fS) 

J B~(S)e v x * m dV 


(29) 


m 


The group effective transport equation for a gas with or without gradients and with a 
cold wall at L is 
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I 


*m = 


B m (S) e 




SmN'e.m}® 


= /' T e,m (L) B e ‘ T e,m d 
J 0 m e, m 


(30) 


The integrated intensity may be obtained by summation of equation (30) over all a 
groups. 

a 

I-E’m (3D 

m=l 


Exact classes of p roblem. - Equation (30) is exact for the following two additional 
classes of problem besides the three classes mentioned subsequent to equation (22). The 
first additional class of problem is the class where the elementary diffusion approxima- 
tion (ref. 1) is applicable and where the wave-number group is narrow enough that the 
variation of B~ between v m and may be neglected. This exactness is proven 

in appendix A where it is also shown that the same radiant flux density is obtained either 
by integrating equation (30) or from the diffusion approximation. In the second additional 
class of problems where equation (30) is exact, a~ must be of the form 


a ?(S) = w m (?)t m (S) 


(32) 


for v between V and ^ m+ p and the wave-number group must again be narrow 
enough that the variation of B~ between V and may be neglected. When these 

two conditions are satisfied, the exactness of equation (30) is proven in appendix B (eq. 
(Bll)) for cases with or without temperature or pressure gradients. 


Refinement by Using Origin Groups 

Error due to a~ not factoring. - The underlying assumption in the application of 

equation (30) to general problems involving temperature and pressure gradients is that 
a~ may be factored approximately into the form given in equation (32). This assumption 
may break down in two types of problems: 
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(1) Problems involving inhomogeneous mixtures 

(2) Problems involving large ratios of maximum to minimum temperature together 

with a gas possessing spectra of widely different excitation potentials 
Introduction of origin groups . - Two spectral absorption coefficients a~ . and 
a~ are said in this report to be uncorrelated within wave -number group m if 


B v^v, j " a Pl, m, j^ a ?, k " a Pl, m, “ 0 


(33) 


The quantities ap^ m j and ap^ m ^ are group Planck mean absorption coefficients 

for a~ . and a~ ., respectively, and are defined by 
v >s 


a Pl, m, v 



B„a, 


V A V,v 


dv 


(34) 


where v can be j or k. 

If spectra are sufficiently complex, then little or no correlation generally exists be- 
tween a~'s for different gases of a mixture or, for that matter, between the low- 
excitation -potential and high-excitation-potential contributions to a~ of a single gas. 
This lack of correlation suggests the approach for refining equation (30). In many cases, 
a spectrum for which equation (32) is not even approximate may be so divided into u 
origin groups r = 1, 2, 3, . . . u that there is little or no correlation between the 
spectral absorption coefficients of any two origin groups, but the spectral absorption 
coefficient a~ r of each origin group approximately factors the way a~ does in equa- 
tion (32). When origin groups are used, the spectral absorption coefficient a~ of the 

gas or mixture is the sum of all the a~ . For an example of the use of origin groups, 

is, r 

if an inhomogeneous mixture contains three gases a, b, and c, and only gas c has spec- 
tra of widely different excitation potentials, then origin group 1 might correspond to all 
the spectrum of gas a, origin group 2 to all the spectrum of gas b, origin group 3 to the 
low -excitation-potential spectrum of gas c, and origin group 4 to the high- excitation- 
potential spectrum of gas c. 

If no correlation is assumed to exist between the a~ of different origin groups, 
then there should be little or no preferential absorption of, say, radiation of origin 
group 1 by origin group 2. Consequently, a g m should be replaced by the sum of the 
a e m r of the various origin groups, and the value of each a g should depend on 

the group effective optical depth r for that Origin group and no other. This de- 

pendence is obtained by defining the effective absorption coefficient of the r origin 

i.u 

group in the nr wave -number group by 
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(35) 


a e, m, rl*®- S i, m, r] * 


/ -^Ti ^) S i rv-» „ 

B v ty)a. v (i//)e ’ ’ ’ dv 

j_ 

/ -a~ _(i^)S. „ _ 

B~(i//)e r x » m ’ r dv 


Here xf/ now stands for temperature and partial pressures of all chemical species at S 

(which, in general, all affect the width of spectral lines of origin group r). The quantity 

S. is defined as the distance along a ray in an isothermal, isobaric gas having only 

l, in, r 

the absorption coefficient a~ and having temperature and pressure the same as at S, 

such that the effective optical depths T e m r of the r^ 1 origin group in the m** 1 wave- 

number group in the isothermal, isobaric gas and in the gas with gradients are the same. 

In this definition both gases are assumed to have the same line widths. 

The group effective optical depth r is given by 

e, m, r 




m, r 


(S')]dS' .T e m r = ln 


B m< S > 


l 


B~(S)e 




(36) 


The total group effective optical depth r g m is 


u PS u 

;, m “ ^ T e, m, r _ / ^ a e, m, r[^^® ^ r e, m, r]^® 

r=l J 0 r=l 


(37) 


The group effective transport equation for a gas with or without gradients and with a 
cold wall at L is no different from equation (30) except for the evaluation of T e 


'-■f 


e » m ^B (t 

m v e,m 


)e" Te » m dr 


e, m 


(38) 


It can be shown that the use of origin groups as outlined is equivalent to the assumption 
that radiation belonging to one origin group is absorbed by another origin group in the 
same proportion as radiation from a blackbody would be. 

Exact classes of problems . - Equations (37) and (38) are either exact or an excellent 
approximation for six classes of problems. 

The first class is a homogeneous mixture of any number of slightly nongray (slightly 
colored), isobaric gases with or without temperature gradients and with a cold, black 
wall behind it. The spectral absorption coefficients of the various gases are assumed to 
be uncorrelated and independent of position. For such a mixture, a~ r of each gas r 
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may be represented by 


a~ = a 

g,m, 


+ Aa 
r r 


(39) 


where a g m r is the average value of a~ r in the wave-number group, and Aa f is any 
small wave -number- dependent perturbation. If the wave-number group is narrow enough 
that the variation of B~ between V m and ^ ^ may be neglected, then both the exact 

group transport equation obtained by integrating equation (3) assuming a cold wall 

/* /*Tpj(L) —7 Vj 

(40> 


and the effective group transport equation (eq. (38)) may be expanded in series involving 
integrals of powers of Aa^a^ m where a g m is the sum of all a g m r with the same 
m value. It is shown in appendix C that the two series expansions so obtained are iden- 
tical through third order in Aa /a , which proves that equations (37) and (38) are an 
excellent approximation for this class of problem. 

The second class is an inhomogeneous, isothermal mixture of any number of slightly 
nongray gases that may possess pressure gradients and have a cold, black wall behind 
them. The spectral absorption coefficients of the various gases are assumed to be un- 
correlated. For such a mixture, a~ r of each gas may be represented by 


= (w„ „ + Aw Jt (S) 

v, r g,m,r r 7 m,r x 1 


(41) 


where w g m r is a constant for the particular gas and wave number group, Aw r is any 
small perturbation depending only on wave number and having an average value of 0 in 
the wave -number group, and t r is a factor depending only on S and due to concen- 
tration variation. Assuming the wave-number group may be any width, both the exact 
and the effective group transport equations may be expanded in a series involving inte- 
grals of powers of Aw /w . It is shown in appendix D that the two series expan- 

sions thus obtained are identical through third order in Aw /w , which shows that 

equations (37) and (38) are an excellent approximation for this class of problem. 

The integrated or group intensity for class 3, consisting of a homogeneous, iso- 
thermal, isobaric mixture of any number u of gases with lines possessing Lorentz con- 
tours and with no continua, may be calculated by the Mayer-Goody statistical model 
(refs. 7 to 9) if the wave numbers of the lines have a random distribution, if the spectra 
of the u gases are uncorrelated, and if the integrated linear absorption coefficients S r 
of all lines of each gas meet two requirements. The integrated linear absorption coeffi- 
cient of a line of gas r is defined by 
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(42) 


■/ 

Jo 


a P", line ^ 


where a~ 1 - n p is the linear absorption coefficient due to only the line being considered. 
The first requirement on S r is that the probability dS r of ~ 
being between S and S + dS is 


S r of a line of gas r 


P r dS 
r r 



(43) 


— 2 
where S r is the average integrated linear absorption coefficient (in units of, say, cm ) 

for all lines of gas r. The second requirement is that there is no correlation between 
the values of S r and the wave numbers of the respective lines. The gases may be 
present in any proportions. K the half -widths of the lines and the mean line spacings of 
all the gases are the same, then the Mayer-Goody statistical model and equations (37) 
and (38) give the same integrated or group intensities for all depths L of the gases, as- 
suming the gases are bounded by a cold, black wall (fig. 1). Hence equations (37) and 
(38) are as good an approximation as the Mayer-Goody statistical model for this class. 

The fourth class of problems may be called "independent, strong line" problems. 

If two or more gases have very narrow spectral lines and no continua and if the depth L 
of the gases is neither too small nor too large, the centers of the lines may be optically 
thick, but the wings of the lines may be optically thin and cover most of the spectrum. 
Because the coincidence of optically thick line centers of different gases is quite rare for 
gases with uncorrelated spectra, the integrated intensity of the mixture should be the 
sum of the integrated intensities due to each gas, or, from equation (30), 


I 

m 



B (t 
m v c 


> 


-T 


e,m, r 


dr 


e,m, r 


(44) 


where the mixture may possess temperature and pressure gradients. However, if the 
wings of the lines are optically thin and cover most of the spectrum, then t & f « 1, 
and equation (44) becomes 


m 



B m dT e,m,r 


( 45 ) 


Equation (45) may be compared to equations (37) and (38) if it is assumed that 
t « 1; thus, equations (37) and (38) give 

v j All 
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(46) 


m 



a 

e,m,r 


dS 


This equation is equivalent to equation (45). Thus, equations (37) and (38) give correct 
values of for independent, strong line problems. 

The last two of the six classes of problems are (5) mixtures of gray gases and (6) 
mixtures of gases such that the mixture is optically thin at all important wave numbers. 
For both classes, the mixtures may be inhomogeneous and have temperature and pres- 
sure gradients. For these classes, equations (37) and (38) give exact values of I m . 


Allowance for a Black Wall 

Extension of approximations . - Extension of the approximations given in the pre- 
ceding three subsections to problems where the wall is black but not cold is not difficult. 
If a wave-number group is narrow enough that the variations of the Planck functions of 
both the gas and the wall can be neglected between and and if a~ is given by 

equation (32), then it is shown in appendix B (eq. (BIO)) that the group intensity is given 
exactly by 


m 


=f 


(L) -T 

e ’ m B (r )e 
m' e,m' 


e,m dr + B e 
e,m m,w 


- T e,m< L > 


( 47 ) 


even in the presence of temperature and pressure gradients. Here B m is the group 
Planck function for the wall temperature. When the spectrum is not divided into origin 
groups, equation (47) may also be written 



+ B e 
m,w 


-f 

J o 


a 


e,m 




(48) 


Substitution of equation (37) into the second term of equation (47) reveals that the use of 
origin groups is equivalent to assuming the multiplication property of transmissions 
(ref. 10). This property has been observed experimentally for molecular band spectra 
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of different gases (ref. 10). 

Exact classes of problems . - Equation (47) is exact for four classes of problems: 

(1) the class previously discussed for which a~ is given by equation (32), (2) gray gases 
or mixtures of gray gases, no matter how wide the wave-number group is made, (3) 
gases or mixtures that are optically thin at all important wave numbers, irrespective of 
how wide the wave-number group is, and (4) the trivial case of an isothermal, isobaric 
gas with a black wall at the same temperature as the gas and with any width of wave- 
number group. 

Hypothesis for other problems . - The hypothesis is made here that equation (47) may 
be applied as an approximation to cases not falling within the four classes of problems 
previously given. Support for this hypothesis is provided by specific cases in the follow- 
ing section. 


RESULTS AND DISCUSSION 
Behavior of the Effective Absorption Coefficient 

The general behavior of the effective absorption coefficient has been discussed in the 
section Behavior of a g for a homogeneous, isothermal, isobaric gas (p. 10). Here a 

specific example, profile A, is considered. Profile A is a spectral absorption coeffi- 
cient consisting of the step function shown in figure 2 and tabulated in table I. The be- 
havior of a g for profile A is shown for an isothermal, isobaric gas in figure 3 together 
with the behaviors of the Planck mean absorption coefficient apj and the Rosseland 
mean absorption coefficient ap Q (ref. 1) defined by 



All three mean absorption coefficients always lie between the minimum and the maximum 
spectral absorption coefficients a^ min and a^ max , respectively. Where the spectral 
absorption coefficient and the pressure and temperature do not vary with S, as in the 
present case, apj and ap Q are independent of S, but a g is dependent on S. The ef- 
fective absorption coefficient a g decreased from apj to a~ m ^ n as S varied from 
zero to infinity. 
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The effect of temperature on a g is shown in figure 4 for an isothermal, isobaric 
gas. The dependence of a g on temperature was caused by the shift of the Planck function 
as shown in figure 2. 

The effective optical depth for an isothermal, isobaric gas with profile A is also 
shown in figure 4. The effective optical depth did not depend linearly on S but increased 
more slowly at high values of S due to preferential self -absorption. 


Methods of Evaluating Integrated Intensity 

Integrated intensity was calculated for specific cases by several approximate meth- 
ods. 

(1) The ’’apj" method was evaluation of the gray-gas group transport equation using 
the group Planck mean absorption coefficient apj m in place of a g m [^/ (S) , T g m J and 
a ^(S T ), T* 1 in equation (48). The equation for a pl is equation (34) without the 
subscripts v. This procedure may be considered a first approximation to equation (48) 
because it gives exact values of I m for gases that are optically thin at all important 
wave numbers. 

(2) The "one-origin-group a " method was the evaluation of equation (48) and thus 
is equivalent to doing a problem without introducing origin groups. 

(3) The "two-origin-group a e " method was the evaluation of equation (47) using 
r _ from equation (37) with u = 2. 

6 y m 

Variations of these methods involved the use of different numbers of wave-number groups. 
For comparison, integrated intensity was also calculated exactly by equation (4). All 
calculations were done on the Lewis IBM 7094 computer. 


Specific Cases with One Gas 

Thirty-six specific isobaric cases involving eight different a~ profiles were se- 
lected to test the accuracy of the one-origin-group a g method. Thirty of these cases 
did not belong to any of the five classes of problems for which the one-origin-group a g 
method is exact. The eight a~ profiles are given in figures 2 and 5 to 8 and are tabu- 
lated in table I. Profiles B and C (fig. 5) are constant half-width Lorentz lines, profile D 
(fig. 6) is a step function, and profiles E and F (fig. 6) are linear variations. Profile G 
(fig. 7) is a step function selected to cause a very large error in I when only one wave- 
number group is used. Profile H (fig. 8) is a Lorentz line whose width varies as T - */^ 
as would be expected in many practical, isobaric problems. The arbitrarily assumed 
temperature distributions are designated 1 to 8 in figure 9. A black wall was located at 
L = 50 centimeters in all cases. 
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The integrated intensities calculated by the exact, one-wave-number-group a pp and 
one-wave-number-group, one-origin-group a g methods are given in tables n and m for 
profiles A to F. The ratios of approximate to exact integrated intensities are plotted in 
figure 10 for profile B. For profiles A to F the a p j method gave errors in I of -43 to 
455 percent, whereas the one-origin-group a g method demonstrated its superiority, at 
least for Lorentz lines, by giving errors in I of only -25 to 24 percent. 

The calculated values of I for profile G are given in table IV. Both the a pl and 
the one-origin-group a g methods with one wave-number group resulted in very large 
errors in I. The difficulty was remedied by using two wave-number groups. Because 
the gas was gray in both groups, the errors vanished for both approximations as shown 
in table IV. 

The values of I for profile H are also given in table IV. Because of the large ratio 
of maximum to minimum temperature, one wave-number group for the one-origin-group 
a g method gave an error of 35 percent. When three wave-number groups were used, this 
error was reduced to 27 percent. The a p j method was totally unsatisfactory with either 
one or three wave-number groups. 


Specific Cases with Gas Mixtures 

Three specific, isobaric cases with mixtures involving five hypothetical gases were 
selected to test the accuracy of the two-origin-group a g method. None of the cases be- 
longed to any of the six classes of problems for which the two-origin-group a g method 
is exact or can be shown analytically to be an excellent approximation. 

The five gases included one completely transparent gas c and four gases a, b, d, 
and e with Lorentz lines. Each of the latter gases was assumed to have 20 lines with 
various integrated absorption coefficients and random wave numbers of the line centers, 
as discussed in appendix E. Consequently, the spectra of gases a, b, d, and e were es- 
sentially uncorrelated. In mixtures, the contribution of gas r (a, b, d, or e) to a~ was 
taken equal to the mole fraction X r times the value of a~ for pure gas r. 

The case referred to as binary mixture A involved gases a and b distributed as 
shown at the right of figure 11. A cold, black wall was located at 50 centimeters, and 
temperature distribution 8 (fig. 9) was chosen. The a~ of the mixture is shown in fig- 
ure 11. The calculated values of I are given in table V. 

The case referred to as binary mixture B involved gases d and e in place of gases a 
and b, respectively, so that the a^, of the mixture (not shown) was different. Other- 
wise binary mixture B was the same as binary mixture A. The calculated values of I 
are given in table V. 

The case named ternary mixture was obtained by inserting a completely transparent 
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gas c between gases a and b of binary mixture A as shown at the right of figure 12. The 
curves in figures 11 and 12 were plotted with the Lewis CACOMP digital computer pro- 
gram and a CALCOMP plotter. Insertion of the transparent gas increased the depth of 
the gas to 150 centimeters, to which point the cold, black wall was relocated. Tempera- 
ture distribution 9 was chosen because it had the same maximum and minimum tempera- 
tures as for the binary mixtures. The a~ of the mixture is shown in figure 12. The 
calculated values of I are given in table V. Despite first impressions conveyed by fig- 
ures 11 and 12, binary mixture A and the ternary mixture had appreciable self-absorption 
(t~ not negligible compared to 1) at most wave numbers. 

A review of the accuracy of the various approximate methods for calculating I (ta- 
ble V) for the three mixtures shows that the a pl method with one or three wave-number 
groups was in some cases unsatisfactory. Calculations by the one-origin-group a g and 
the three -origin- group a g methods were performed with three wave -number groups. 

In these calculations, one origin group gave errors between -49 and 19 percent, while 
two origin groups gave errors between -17 and 28 percent, demonstrating the superiority 
of the two- origin- group a g method for inhomogeneous mixtures containing two semi- 
opaque gases. 


RECOMMENDATIONS FOR PRACTICAL APPLICATIONS 

The approximations developed in this report are well suited to calculations with a 
high-speed digital computer but not well suited to hand calculations. Because the effec- 
tive absorption coefficient approximation introduces errors of the order of ±28 percent in 
the calculation of integrated intensity and hence radiant flux density, the numerical meth- 
ods used in practical calculations can be crude without introducing serious additional 
error. The proposed procedure for practical calculations of radiant flux density is out- 
lined in the following sections. 


Selection of Wave-Number Groups 

Once the gases and their spectra involved in a given problem or series of problems 
have been determined, the next step is to select the wave-number groups. Because 
equation (47) is exact for gray gases, errors in I will be minimized if the spectrum is 
so divided that each wave-number group is as nearly gray as possible. A very conserva- 
tive approach would be to divide the spectrum so that no Planck function corresponding to 
any temperature of the problem varied appreciably in any wave-number group. When a 
large number of closely related cases are to be calculated, a better approach would be 
to calculate a few cases exactly by applying equation (4) to the entire spectrum and then 
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to determine by trial and error how many wave-number groups would be needed if equa- 
tions (31) and (47) are used as an approximation. 


Selection of Origin Groups 

For problems involving chemical reactions, at least one origin group generally will 
be required for each chemical species possessing a spectrum in a given wave-number 
group. 


Calculation of a e m r and x e m r 


The effective absorption coefficients a g m r can be calculated once and for all for 
a series of problems involving the same chemical species. The computer can then store 
them as multivariant polynomial fits or in multiargument tables. If one origin group is 
assumed per species, it will be necessary to have at most u + 2 arguments (T, m r , 
and partial pressures of all u chemical species). However, in many cases fewer argu- 
ments will suffice because the various particles that cause significant collision broaden- 
ing of spectral lines often produce the same line shape and because a g m r /p r may be 
stored instead of a g m r itself (p r is the partial pressure of species r). Thus, 

a /p often may be stored with the three arguments T, r , and b where b 

e,ni,r r ^ , m , x u c 

is an average collision half-width. 

It was mentioned in the section Nongray Gases with Gradients that in constructing a 
table of a g values, a number of values of S i are substituted into equations (17) and (19) 
to get a and r . Similarly, in constructing a table of a values, a number of 

values of m r may be substituted into equation (35) and the right-hand side of equa- 
tion (36). A better procedure is to nondimensionalize equations (35) and (36) by defining 
dimensionless photon energies by 


y = he? 
kT 
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hcV 
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kT 
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(50) 


From equations (2), (28), 


(36), and (50), 
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(51) 


e, m, r 


= In 



dV 



dV 


where a^, r has been rewritten ay r for clarity, and m under the integral sign means 
to integrate from V m to V m + r Evaluation of the integral in the numerator can be 
rapidly carried out using the expressions for Q(V) given in appendix F. From equations 
(2), (35), and (50), 


e,m, r 



(52) 


For spectra involving many lines, use of statistical models (refs. 7 to 9) may be 
preferable to use of equations (51) and (52). 


Evaluation of Integrated Intensity 

If one origin group is assumed per species, integrated intensity I m may be evalu- 
ated at point P for a ray in a given direction (fig. 1) in an inhomogeneous mixture of u 
gases by integrating by numerical methods u first-order differential equations and a 
single integral. The u first-order differential equations are given by 


dr. 


e>m,r = a 

dS e ’ m 


, r[^^ ,T e,m,r] 


(r = 1, 2, 3, . . ., u) 


(53) 


which is equivalent to the left-hand side of equation (36). The single integral is given in 
equation (47), in which B m may be evaluated by the method in appendix F. 
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Evaluation of Radiant Heat -Flux Density 


The component F n of the heat-flux density at the point P in the direction n nor- 
mal to the infinitesimal area dA (fig. 1) may be calculated from the relation 

/a 

£ I m cos 0 d a (54) 

m=l 

where dS2 is an element of solid angle with the ray at its center, 6 is the angle between 
the ray and the direction n, and the integration is performed over 4 n ste radians. 


CONCLUSIONS 

1. The effective absorption coefficient approximation simplifies radiant-heat trans- 
port calculations for cases with complicated spectra that are optically thick at some im- 
portant wave numbers. Integration over all important wave numbers is performed once 
and for all for all gas mixtures containing the same chemical species. 

2. To obtain accurate integrated intensities where widely different temperatures oc- 
cur in a problem, it may be necessary to divide the spectrum into wave-number groups, 
with an effective absorption coefficient and an effective transport equation for each 
group. 

3. The use of origin groups in the calculation of the effective absorption coefficient 
is generally necessary to obtain accurate integrated intensities for inhomogeneous mix- 
tures and may be desirable where spectra of the same gas have widely different excita- 
tion potentials. For valid results, spectra of the different origin groups must have no 
appreciable correlation with each other, and the spectral absorption coefficient of each 
origin group must factor approximately into wave-number-dependent and position- 
dependent functions . 

4. When one to three wave-number groups were used as required and one or two 
origin groups were used as appropriate, the effective absorption coefficient approxima- 
tion gave integrated intensities with errors of -25 to 28 percent for 39 cases. In con- 
trast, use of Planck mean absorption coefficients in the gray-gas transport equations for 
the same number of wave-number groups gave integrated intensities with errors of -43 
to 455 percent. 
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5. If hot, gray or black radiation falls on a cold gas with narrow, isolated spectral 
lines, the fraction of the radiation transmitted may be underestimated seriously by using 
the Planck mean absorption coefficient. The error can be reduced substantially by using 
the effective absorption coefficient approximation. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, February 16, 1967, 

122-28-02-17-22. 
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APPENDIX A 


DERIVATION OF GROUP DIFFUSION EQUATION FROM 
EFFECTIVE GROUP TRANSPORT EQUATION 

It is shown in this appendix that if the wave -number group is narrow enough and if the 
usual assumptions of the elementary diffusion approximation (ref. 1) are made, then the 
effective group transport equation gives exact values of I m< The group diffusion equation 
is derived from the equation for I m and is found to be the same as derived from the 
literature. 


EXACTNESS OF EFFECTIVE GROUP TRANSPORT EQUATION 


Consider a gas where the Planck function does not vary appreciably in a few photon 
mean free paths at any wave number. Then the elementary diffusion approximation 
(ref. 1) should, in general, be applicable. Although it is not obvious from the usual der- 
ivation (ref. 1), this approximation entails assuming that a~ is independent of tempera- 
ture and position. This assumption is reasonable because when the approximation is ap- 
plicable all significant contributions to the radiant flux density come from gas within a 
few photon mean free paths of the point in question, and the Planck function (and there- 
fore temperature) must not vary appreciably in this region. Hence, in the following der- 
ivation, is assumed to be independent of temperature and position. 

If a wave-number group is chosen narrow enough that the variation of. B~ with V 
can be neglected, then according to equation (27), a g m is given by 
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Substitution of equation (Al) into equation (29) and differentiation of the result with re- 
spect to S give dS. /dS = 1. By using the boundary condition S. = 0 at S = 0, it 

X y 111 1 J III 

is found that S i m = S for all values of S, and equation (Al) becomes 
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Thus, 
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If the gas is assumed to be infinite, equation (30) becomes 


(A3) 


(A4) 



where q m (0) is the value of q m for S = 0, and is 1/ - V . From equa- 

tion (A3), 


/ _ a ^s 
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Combining equations (A5) and (A6) gives 
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(A6) 


(A7) 


which is exact because it may also be obtained from equation (3) by setting L = °°, 
dropping the last term, integrating over the wave-number group, and neglecting the 
variation of B~ between V m and . 
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DERIVATION OF GROUP DIFFUSION EQUATION 


The group diffusion equation may be obtained from equation (A7) if B m is expanded 
in a Taylor series about S = 0. 

j + . . . (A8) 

C 

and only the first two terms are retained. Combining equations (A7) and (A8) and 
changing the order of integration give 
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'tiB 

= B (O) +S(— “ 
m 1 dS 



/dB„ 

B - ( ° )tS br 


O 


e"^ 

dr~ diX 


AV 


(A9) 


m 


This equation simplifies to 


/dB m 
I = B (O) +( — - 
m m ' 7 \ dS 


I A-!--* 

O / a~ AV m 
m 


(A10) 


From Stewart and Pyatt (ref. 11), the local Rosseland mean absorption coefficient 

a— for the wave-number interval AT? is defined by 
Ro, m m 


— i — * f-L- 

a Ro,m J m a ? A "m 


(All) 


where the variation of dB~/dT in this wave-number interval has been neglected. Com- 
bining equations (A10) and (All) gives 


a Ro,m\ /O 


(A12) 


The gradient of T is assumed to be in the direction n or -n. The group heat flux 
density is (fig. 1) 


m 


= - / I cos 0 dfin 

*4r m 


(A13) 
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where, from equation (A12), 


I = B (O) + 
in m 


cos 6 
a Ro,m 



(A14) 


Here x is the component of S in the direction n. Combining equations (A13) and (A14) 
gives 



Ro, m 


(A15) 


This is simply the group diffusion equation, which can also be derived from the spectral 
diffusion equation (ref. 1) by integrating it over the wave-number interval with the 

assumption that (dB~/dx)Q is constant in this interval. 
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APPENDIX B 


DERIVATION OF EFFECTIVE GROUP TRANSPORT EQUATION FROM EXACT 
GROUP TRANSPORT EQUATION FOR A CLASS OF PROBLEMS 

In this appendix the effective group transport equation is derived from the exact 
group transport equation for a gas whose a^, can be separated into wave-number- 
dependent and position- dependent factors (eq. (32)). Temperature and pressure gradients 
may be present. The gas is bounded by a black wall at S = L, which may be at any tem- 
perature. The wave -number group must be narrow enough that the variation of the Planck 
functions of the gas and of the wall for wave numbers between V and V ^ may be 
neglected. 

Equation (26) is the exact group transport equation. If the variation of B~ and 
B~ w between V and ^ is neglected and the order of integration is changed, 
equation (26) becomes 


1 /* L f -T w( s ) B f -t w (L) 

1,^=— / B (S) J e V z v (S)dVdS + J^L J e V dV (Bl) 

AV m AV m Jm 


The next step is to find a simple expression for the product M S)S i m( S >- The 
quantity m (S) is defined by equation (29). 

If the variation of B~ between V m and ^ m+ j is neglected, equation (29) becomes 


X Sa e,m[ S '> S i,m( S ')] ds '= I ” 


AP! 


m 


/ 

•'m 


-a~(S)S. (S) 

x v v ' i,m v ’ 


dv 


(B2) 


Equation (27) defines a. . If the variation of B w between PL and V 1 is ne- 

e, m v m m+1 

glected, equation (27) becomes 


a e,m[ S ’ S i,m< s >] * 


f , x -a«(S)S (S) 

/ a~(S)e v 1 » m dV 
y m v 


(B3) 


Equation (B3) is then substituted into equation (B2), equation (32) is substituted into the 
result, and the equation thus obtained is differentiated with respect to S. 


34 



dS. IS) S. (S) dt (S) 
i,m x ' i,m x _ m v 7 _ ^ 

dS 


t m (S) dS 


(B4) 


The solution of the differential equation (B4) with the boundary condition = 0 at 

S = 0 is 


/ 

Jo 


/a m 
S. (S) =— 

1 ’ m t ( 

m v 


t m (S’)dS’ 


(S) 


(B5) 


Combining equations (32) and (B5) gives 


V(S) s i>m (S) = a ~(S , )dS* = 7y 


(S) 


(B6) 


This equation gives the desired simple expression for a~(S)Sj m (S). 

Substituting equation (B6) into equation (Bl) gives 

< /*L r -a~(S)S. m (S) B w r -a~(L)S. m (L) 

I = — - — / B (S) / e v 1 » m a~(S)dy dS + m ’ w / e v 1,m 

m Av_ Jo m *^n v AV 

— m 


d? 


(B7) 


If the variation of B^ between and is neglected, the right-hand side of 

equation (29) gives 


e - T e, »«__!_ /■.-MW 1 . m « dr 
AS, m m 

If equation (B8) is multiplied by equation (B3), the result is 


(B8) 


e -T e,m l ® ) (s) =t — — / V 

e,m Atr * / m 
m 

Equations (B8) and (B9) are then substituted into equation (B7) with the result 


(B9) 
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/ ^ -T e m® 
B m (S)e e ’ m 


E (S)dS + B g 
e,m v ' m,w 


■ T e,m< L > 


(BIO) 


The derivation of equation (BIO) has shown that it is an exact equation for the class of 
problem considered. If the wall is cold (B negligible), equation (BIO) becomes 

JI1 2 w 

/ L -r (S) 

B m (S)e a em (S)dS (Bll) 

Equations (BIO) and (Bll) are effective group transport equations. 
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APPENDIX C 


SERIES EXPANSIONS OF TWO GROUP TRANSPORT EQUATIONS FOR 
HOMOGENEOUS MIXTURE OF NONGRAY GASES WITH OR 
WITHOUT TEMPERATURE GRADIENTS 

In this appendix, series expansions of the exact group transport equation and of the 
effective group transport equation (eq. (38) with T e>m given by eq. (37)) are carried 
out and compared in order to help establish the validity of the effective group transport 
equation with more than one origin group. 

The class of problem considered is restricted to a homogeneous mixture of u 
slightly nongray (slightly colored) isobaric gases in any proportions and with or without 
temperature gradients. The gases are assumed to have uncorrelated spectra and are 
numbered 1,2, . . . , u. The two expansions give the group intensity at S = O (point P 
in fig. 1). The distance L along the ray to a cold, black wall can be any value including 
infinity. The optical depth of the gas along the full length of the ray is arbitrary. The 
wave-number group considered is assumed to be narrow enough that the variation of the 
Planck function B~ with wave number can be neglected for the range of wave numbers 
included in the group. 

The spectral absorption coefficient of each gas of the in wave-number group is 
assumed to be independent of position along the ray and may be expressed as 


= a 


g,m,r 


+ Aa 


r = 1, 2, . , u 


(Cl) 


where r designates the particular gas. The gray absorption coefficient a is a 

constant for each wave-number group for each gas, so the perturbation Aa r is some 
function of wave number. Integrating equation (Cl) over the range of wave numbers in- 
cluded in the wave-number group m gives 




r 


dV + 



Aa df 
r 


or 

f m a *, r dlr “ a g, m, r * j [ “r di? 
If a^, r is given, a g m r may be so selected that 


(C2) 
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diX = a ASX r = 1, 2, . . . , u 

r g,m,r m 7 7 7 


f a ?, 

' / m * 

Subtracting equation (C3) from equation (C2) gives 


(C3) 


l 


Aa dV= 0 r = 1, 2, . . . , u 


(C4) 


The requirement that each gas be slightly nongray may be specified by 

Aa, 

« 1 


(C5) 


I a g,m,r| 

for all wave numbers in the wave-number group. The total spectral absorption coeffi- 
cient is given by 


r=l 

Equations (Cl) to (C6) apply in both the series expansions which follow. 


EXPANSION OF THE EXACT GROUP TRANSPORT EQUATION 

From equations (1) and (26), the exact group transport equation with a cold, black 
wall is 


m 






dS diX 


(C7) 


Because it was assumed that the variation of B~ with wave number could be neglected, 
equation (C7) may be rewritten 


m 



" a ix s 


a^diXdS 


(C8) 


If equations (Cl) and (C6) are used, equation (C8) becomes 
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Define a total gray absorption coefficient for wave-number group m by 


u 


a g,m 2 a g,m,r 
r=l 


Also, define a gray optical depth by 


(CIO) 


t = a S 
g* m g,m 


(Cll) 


Combining equations (C9), (CIO), and (Cll) and expanding the Aa^ terms of the exponen- 
tial in powers of Aa^/a^ m give 



where r (L) is the value of r at S = L, and TT indicates a product. 

g,m' 7 g,m >11 

Appropriate expressions will now be found to express the assumption that the spec- 
tra of the u gases are uncorrelated. The perturbation Aa r at a given wave number 
may be associated with a hypothetical photon of the same wave number, the photon being 
a member of a hypothetical population r of photons with a constant frequency function 
between V and and a zero frequency function outside these limits. In other 

words, the probability of a photon having a wave number in an infinitesimal range dV 
'of fixed width is independent of wave number between and *m + l and is zero out- 
side this range. K there are u statistically independent populations of photons with this 
type of frequency function, then the following averages exist (ref. 12). 
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dV= 0 


j*kl 



Aa^ Aa k 


Aa^ = 0 


j*k 
j * l 

k * l 



3 * k J 


(Cl 3) 


Equation (C13) defines "u gases with uncorrelated spectra. " It is a good assumption for 
atomic and ionic spectra consisting of many lines by themselves or superimposed on con- 
tinua which do not have strong wave-number dependence. 

Equations (C4) and (Cl 3) reveal that many of the terms that occur in evaluating 
products in (C12) are zero, with the result that 




(Cl 4) 


where terms higher than third order in Aa /a have been dropped as insignificant. 
Equation (C14) is the truncated series expansion for the group intensity at point P based 
on the exact group transport equation. 


EXPANSION OF THE EFFECTIVE GROUP TRANSPORT EQUATION 

From equations (37) and (38) ? the effective group transport equation may be written 


m 



-T 


e, m, r 


dr 


e,m 


(Cl 5) 
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Consequently, evaluation of equation (Cl 5) requires two preliminary steps: the evaluation 
of d'r (from which dT is obtained) and the evaluation of exp(-r ). 

G y IH ) 1 0 , HI " ) ill J * 

The value of a Q _ is needed to evaluate dr .If the variation of B~ with 

G y m y r G y III j 1 V 

wave number within the wave-number group is neglected, equation (35) becomes 


f *v,. 

*'m ’ 


'V, r S i, m, r 


e,m, r 


/-a— S . 

e u > r 1,m,r dP r 


m 


(Cl 6) 


Substitution of < lation (C16) into equation (36) and differentiation of the result with re- 
spect to S giv< dS. /dS = 1. If the boundary condition S. = 0 at S = 0 is 
used, it is found that S. = S for all values of S. From equations (C16) and (36), 


dr 


* / m 


a~ r dVdS 


e ’ m , r r -a~ S 

/ e u > r dV 
*"m 


(Cl 7) 


Let the denominator of equation (C17) be designated by 


Vr 



diX 


(CIS) 


Substituting equation (Cl) into equation (C18) gives 


%i, r 



S - Aa S 
r r 


d V 


(Cl 9) 


A gray optical depth for gas r may be defined by 

t = a S 

g, m, r g,m,r 


(C20) 


Substituting equation (C20) into equation (Cl 9) and expanding the second exponential in a 
power series give 
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Vr =e 


g,m,r 


A ~ _g 1 m 1 r 

m 2 / \a 


b 


Aa 

— — i dp- g> m > r 

6 


g,m,r. 


/< 


Aa. 


dv + 


a 


g,m,r/ 


(C21) 


where the second term in the series has been omitted because it is zero according to 
equation (C4). 

Similarly, the integral in the numerator of equation (C17) is designated by 


/ -a~ S 
v, r 

e ’ a~ 
v, 

' m ’ 


dV 


(C22) 


Substituting equation (Cl) into equation (C22) gives 


. -f 

m , r J 
7 m 


-SL„ rvi J 3 -Aa S 

g,m ’ r e r (a 


+ Aa )d.p 
g,m, r r' 


(C23) 


Expanding the second exponential in a power series and evaluating the result give 


z = e 
m, r 


g,m,r 


AV + [ S> m > r - t rl 
m \ g g> m > r J 


/ 


Aa. 


g,m,r> 


d? 


„ 2 3 

T T,m,r T g, m, r 1 




Aa. 


g,m,r/ 


dir + . 


a (C24) 

g,m,r 1 ^ 


where the second term in the series has been omitted because it is zero according to 
equation (C4). 

Combining equations (C17), (C21), and (C24) gives 


dT 




dv 


g>m,r 


dS (C25) 
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where terms of higher than third order in Aa r /ag m r have been neglected. 

A different form of dr is useful to evaluate exp(-r ). From equation 

V y III y X ‘Gy 111 y 1 

(Cl 7) and the right-hand side of equation (36), 


dr 


e,m, r 


/ -a~ S 

e v * r a~ „ dFdS 

v,r 


AVe e ’ m ’ r 
m 


(C26) 


This equation can also be written 


e‘ Te ’ m ’ r dT 


z_dS 
_ m, r 


e ,m,r 


m 


(C27) 


Combining equations (C24) and (C27) and trimeating the series at third order in 
Aa r/ a g,m,r S ive 


e" Te ’ m ’ r dr 


-T 


e, m, r 


= e 


g.m.r 


1 + 




g>m,r 



dr 


(C28) 


Integrating both sides from the origin of the ray out to S and simplifying the resulting 
expression give 


e ” T e,m, r _ e _T g,m,r 


1 | g> m > r 


2 AV 


m 



(C29) 

It is now possible to evaluate equation (C15). First equation (C25) is substituted into 
the differential form of equation (37) with the result 
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II 



(C30) 

Substituting equations (C29) and (C30) into equation (Cl 5) and using equations (CIO) and 
(Cll) give 



Evaluating equation (C31) and retaining only third- or lower-order terms in Aa /a 

r ^ j m 

gives equation (C14). Hence, through third order in Aa /a , both the exact group 

r g ? m. ^ 

transport equation and the effective group transport equation give the same value for I m . 
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APPENDIX D 

SERIES EXPANSIONS OF TWO GROUP TRANSPORT EQUATIONS FOR 
INHOMOGENEOUS, ISOTHERMAL MIXTURE OF NONGRAY GASES 

In this appendix, just as in appendix C, series expansions of the exact group trans- 
port equation and of the effective group transport equation (eq. (38) with t q m given by 
eq. (37)) are carried out and compared in order to help establish the validity of the effec- 
tive group transport equation with more than one origin group. 

The class of problems considered is different than in appendix C. It is restricted to 
an isothermal mixture of u slightly nongray gases in any proportions and mixed in a 
homogeneous or inhomogeneous manner. Pressure gradients may be present. The gases 
are assumed to have uncorrelated spectra and are numbered 1, 2, . . ., u. The two 
expansions give the group intensity at S = O (point P in fig. 1). The distance L along 
the ray to a cold, black wall can be any value. The optical depth of the gas along the full 
length of the ray is arbitrary. The wave-number group considered may be any width. 

It is assumed that the absorption coefficient of each gas in wave-number group m 
may be expressed as 


a~ _ = w (tx)t (S) r = 1, 2, 
v , r m,r x ' m,r' ' 7 7 


u 


(Dl) 


where t r is a factor dependent only on S and due to concentration variation. In 
equation (Dl), 


w = w + Aw 

m, r g, m, r r 


(D2) 


where w is a constant, and Aw is a small perturbation, dependent only on wave 

g, m, r ? r 

number. Multiplying equation (D2) by B~ and integrating over the wave-number group 
give 



B^w 


m, r 


dv = B _ _ 
m g,m, r 


+ 



Aw r di? 


(D3) 


If w m r is given, w^ m r may be so selected that 




r d^ = B m w g,m, r r = 2, . . . , u 


(D4) 
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Subtracting equation (D4) from equation (D3) gives 



B~ Aw r dV = 


0 


r = 1, 2, 


u 


The condition that each gas be slightly nongray is specified by 


Aw r 

w 

g,m,r 


«1 


(D5) 


(D6) 


for all wave numbers in wave -number group m. Equations (Dl) to (D6) apply in both 
the series expansions which follow. 


EXPANSION OF THE EXACT GROUP TRANSPORT EQUATION 

From equation (26), the exact group transport equation for an isothermal gas with a 
cold, black wall is 


m 


/ -T~ 

B v 1 " e V 

XXX 


- t w ( l )1 


dF 


(D7) 


From equations (1), (Dl), and (D2), 


A /*L A /*L 

- > w „ / t _(S)dS - > Aw / t 
Zj g»m,ry n m,r' r./ n m, 

r=l u _ r=l u 


-Tp'(L) 
e = e 


Define a gray optical depth m for S = L by 


(S)dS 


(D8) 


A r h 

t (L) s y w / t (S)dS 

g, m' ' /_, g,m,ry n m, r v ' 


(D9) 


r=l 


Substituting equation (D9) into equation (D8) and expanding the second exponential give 
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V L) , »" Tg ’ m(L) ftf - “r/ L ‘m,r (S)dS [** r £ Vr®*] 


j[ iw r/ L ‘ m> rH 3+ -] 


(DIO) 


Substituting equation (DIO) into equation (D7) gives 


= / 




• t s,»i (l) '[T, 

r=l 


1 - Aw 


r 


Jq t m,r^ dS + 2 t m,r^ dS ] 


/•L I 3 

AW r / 


+ . . . 


V dir (Dll) 


In a manner similar to that for equation (C13), the lack of correlation of the spectra 
of the u gases is specified by 


l 


AWj Aw k dV = 0 




j *k 


j * k 


B~ AWj Aw k AWj d"? = 0 ] * l 


l 

K ¥= 6 

y B^ AWj^Aw^ 2 d^ = 0 j * k 


(D12) 


Define a gray optical depth for gas r for S = L by 


g 


/•L 

_ „(L) ■ w (r l t (S)dS 

,m,r' g,m,ry Q m,r v 


(D13) 


With the help of equations (D5), (D12), and (D13), equation (Dll) becomes 
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(D14) 


This equation is the desired expansion of the exact group transport equation. 


EXPANSION OF THE EFFECTIVE GROUP TRANSPORT EQUATION 

From equations (38) and (37), the effective group transport equation for an iso- 
thermal gas is 


I = B 
m m 


■ -T e m (L)‘ 
1 - e e ’ m 


= B. 


m 


u 


" T o n, r.( L )l 


1 -JJ e ' e ’ m ’ r ' 


L r=l 


From equations (35) and (36), 


dr 


i 


B~e" a ^ rSi ’ m ’ r a~ r d^dS 


e,m, r 


B e _T e» m 5 r 


By a derivation analogous to the derivation of equation (B5), 




(D15) 


(D16) 


(D17) 


for all values of S. 

From equations (Dl), (D16), and (D17), 
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dV dS 


-r 


e,m, r 


dT. 



t (S’)dS' 
m, r w 


(D18) 


Integrating equation (D18) from S = O to L gives 


- T 


e,m, r 


(L) 




t (S)dS 
m, r v / 


dy 


(D19) 


Substituting this equation into equation (D15), making use of equation (D2), expanding into 
a power series, and making use of equations (D5), (D9), and (D13) give equation (D14) 
through third order. Hence, through third order in Aw /w , both the exact group 

transport equation and the effective group transport equation give the same value for I . 
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APPENDIX E 


SPECTRAL ABSORPTION COEFFICIENTS FOR GASES IN GAS MIXTURES 

The method of obtaining the spectral absorption coefficients for hypothetical gases 
a, b, d, and e in the gas mixtures is given herein as well as the spectral absorption co- 
efficients obtained. Each gas was assumed to have 20 spectral lines whose wave numbers 
and integrated intensities were given by the Mayer-Goody statistical model (refs. 7 to 9). 


WAVE NUMBERS OF LINE CENTERS 

For each gas the 20 lines were restricted to wave numbers between 0 and 400 000 re- 
ciprocal centimeters because the Planck function for the maximum assumed temperature 
(40 000° K) was negligible at larger wave numbers. The wave numbers V of the line 
centers were assumed to have random values as required by the Mayer-Goody statistical 
model. These values were selected by taking six-digit random numbers from the table in 
reference 13 and discarding numbers higher than 400 000. Rearranged in order of in- 
creasing wave number, these values are given in tables VI and VII for gases a, b, d, 
and e. 


INTEGRATED ABSORPTION COEFFICIENTS OF LINES 

The probability P dS r of the integrated linear absorption coefficient S r of a line 
of pure gas r being between S r and S r + dS r was assumed to be given by equation (43). 
Let f(§_) be the cumulative distribution of S r (ref. 14). 




The inverse function 



is defined by 



(El) 


(E2) 


Hence , 


f" 1 ^) = -S r ln(l - y) 


(E3) 


where y is a dummy variable. Let R be a member of a set of random numbers with 
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rectangular distribution between 0 and 1 and no numbers outside this range. Then values 
of S r with the probability distribution given by equation (43) may be found from the rela- 
tion (ref. 14) 


S r = f -1 (R) = -S p ln(l - R) (E4) 

Let J = §*/S . Then, equation (E4) becomes 

J = -ln(l - R) (E5) 

Because of the rectangular distribution of R between 0 and 1, 1 - R can be replaced 
with R; thus equation (E5) becomes 


J = -In R (E6) 

Values of J for the 80 lines were calculated from equation (E6) by taking 80 six-digit 
random numbers for R from the random-number table of reference 13 and are listed in 
tables VI and VII. 

Two other quantities were needed to specify the Lorentz lines in a pure gas. The 
collision half-width (ref. 7) was chosen as 750 reciprocal centimeters for all lines of all 
gases. The integrated linear absorption coefficient S r , which is roughly an average 
value of S r for all lines of gas r, was chosen as 838 reciprocal centimeters squared 
for each pure gas. 

SPECTRAL ABSORPTION COEFFICIENTS FOR PURE GASES 

The spectral absorption coefficients for pure gases a, b, d, and e were calculated 
using the results from the preceding two sections. The spectral absorption coefficient 
a^ a for pure gas a is the curve at S = O in figure 11, a^. ^ is the curve at S = 50, 
while a.p ^ and a^. are not shown but have the same general appearance as a^,. a and 

a F,b‘ 

Because of the manner in which they were obtained, a~ a , a~ ^ a~ and a~ are 
essentially uncorrelated. 
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APPENDIX F 


GROUP PLANCK-FUNCTION COMPUTER SUBPROGRAM 

Near the end of this appendix is a listing of the group Planck-function subprogram 
written in FORTRAN IV for the IBM 7094 digital computer. This subprogram was used 
in several of the programs used to calculate integrated intensity by different approxima- 
tions and should be useful to other investigators for numerous applications. 

The group Planck function B m was defined by equation (28). Let 

% dsy 

Q(P'jT) =— (FI) 

B 


Then, 

4 

B m (T) = — [ Q ^m’ T > - Q^m+1 ’ T >] (F2) 

7 T 

where a is the Stefan- Boltzmann constant. It is readily shown that 


q<v)=m f “ - y ±- dv’ 

■Vv e V - 1 


(F3) 


where V is given by equation (50). Hence, 

B m (T) = 1.805X10" 12 T 4 [Q(V m ) - Q(V m+1 )] 


(F4) 


Expressions for Q(V) have been given by Wiebelt (ref. 15) and are repeated. Let 
W = V 2 . 


Q(V) = 1 - VW 


0.051330 - v(o. 019249 - vjo. 0025665 - w[3. 0554X10" 5 


- W(5. 6581X1 0 -7 - 1.1 573X1 0' 8 W)]]j 


V < 2 (F5) 
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Q(V) = 0.15399 


5 

Z - {[(MV + 3)MV + 6]MV + 6} V > 2 

M 4 1 J 

M=1 


(F6) 


The subprogram for the group Planck function, based on equations (F4) to (F6), with 
the FORTRAN variable BLACK representing the group Planck function B m follows: 


FUNCTION BLACK ( T ) 

COMMON /B 3/QMEGA 
DIMENSION nMEGA(2>,QUAD<2) 

FAC TOR = 1 • 438 8/ T 
DO 7 1=1, 2 

v=factor*omegaii ) 

IFIV-2.H, 10,10 
10 IF< V-60. ) 2, 11, 11 
l 1 QUAD! I ) = 0 . 

GO TO 7 

1 w=V**2 

QUAD! I )=1 .-V*W*( .05133 0-V* ( .01 9?49-V*{ .0025665-W* ( 3.05 54E-**-W* 

1 ( 5.6581E-7-1. 1573E-8*WH ))) 

GO TO 7 

2 ZEXP=EXP( -VI 
SUM=0 • 

WEXP= 1 . 

no 4 M = l, 5 

TF( V*FL0AT(M)-60. >5,4,4 
5 WEXP=WEXP*ZEXP 
PROD=FLOAT(M ) * V 

SUM = SUM+WEXP*( ( (PR 03* 3. } *PR0D+6. )*PR0D+6. > /FLOAT < M) **4 
4 CONTINUE 

QLI AfM I >=. 1 53 99* SUM 
7 CONTINUE 

BLACK = l • 805E— 12*T** 4* ( QUAD ( l )— QUAD! ? ) 1 

RETURN 

END 

The arguments are actually T, OMEGA (1), and OMEGA (2), although the last two are 
in COMMON. OMEGA (1) is the wave number at the lower limit of the group, and 
OMEGA (2) is the wave number at the upper limit. 

The program is accurate to within 0. 00002 B. 
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TABLE I. - EQUATIONS OF SPECTRAL ABSORPTION 
COEFFICIENTS DESIGNATED A TO H 


Profile a 

Spectral absorption coefficient, 

Range of wave number V for 



which expression is valid, 


cm" 1 

cm" 1 

A 

0. 10 

0 < V< 60 830 


.02 

60 830 <V<^ 

B 

3 767 000 

0<V<oo 


(V - 60 830) 2 + 18 840 000 


C 

376 700 

0 <v < 00 


( V - 60 830) 2 + 188 400 


D 

0.02 

0 < V< 60 830 


. 10 

60 830 < V< « 

E 

6.574X10" 7 v + 0.0052 

0 < v < 00 

F 

-6.574X10" 7 V+ 0. 1048 

0 < v< 159 400 


0 

159 400 < V< « 

G 

2.00 

0 < V < 20 000 


.02 

20 000 <V<™ 

H 

636.6 b 

c 

(ST - 60 830) 2 + b 2 
where b c = 90 000/T*/ 2 

0 < v < <» 


a Profiles are given in figs. 2 and 5 to 8. 



TABLE H- - INTEGRATED INTENSITIES CALCULATED BY THREE METHODS FOR PROFILES A TO F WITH A COLD WALL a 


Profile C 


Exact 

One wave -number group* 3 ap^ 
One wave-number group, ^ one 
origin group a g 


Profile D 


Exact 

One wave -number group* 3 ap^ 
One wave-number group,* 3 one 
origin group a g 


Exact 

One wave -number group* 3 
One wave -number group,* 3 
origin group a g 


a Pl 

one 


Exact 

One wave -number group* 3 a™ 
b 

One wave -number group, one 
origin group a 


a Cold, black wall is at distance L along ray of 50 cm. 
* 3 Wave -number group extended 0 to 400 000 cm'*. 


2.23X10 5 



1. 18X10 5 



6.28X10 4 

1.24X10 6 

455 

5.62 

378 

2.52X10 5 

1.96X10 5 

-12 

1. 18 

0 

6.26X10 4 


Method of calculation 

L 


Temperature distribution 




l 

4 


7 



Integrated intensity, 
L 

o 

W/(cm )(steradian) 

Error from 
exact I, 
percent 

Integrated intensity, 
I, 

W/(cm 2 )(steradian) 

Error from 
exact I, 
percent 

Integrated intensity, 

I, 

W/ (cm 2 )(steradian) 

Error from 
exact I, 
percent 



Profile A 




Exact 

One wave -number group* 3 a^. 

b 

One wave -number group, one 
origin group a g 

1.29X10 6 
1.32 
1. 15 

2 

-11 

5.73X10 5 

6.70 

5.73 

17 

0 

3.30X10 5 

3.90 

3.54 

18 

7 



Profile B 




Exact 

One wave-number group* 3 ap^ 
One wave-number group,* 3 one 
origin group a g 

5.91X10 5 

1.21X10 6 

5.09X10 5 

105 

-14 

3.03X10 5 

5.43 

3.03 

79 

0 

1.59X10 5 
2. 40 
1.57 

51 

-1 


301 

0 


1.08X10 6 



5.73X10 5 



3. 15X10 5 

___ 

1.36 

27 

6.70 

17 

3.46 

10 

1.29 

19 

5.73 

0 

2.99 

-5 


Prc 

file E 




1.21X10 6 

— 

6.05X10 5 



3.08X10 5 



1.45 

20 

6.44 

6 

3. 11 

1 

1.39 

15 

6.06 

0 

2.96 

-4 


Pro 

file F 




1.26X10 6 

— 

6. 35X10 5 

— 

3.65X10 5 



1.32X10 6 

5 

6.72 

6 

3.87 

6 

9.93X10 5 

-21 

6.36 

0 

3.79 

4 
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TABLE m. - INTEGRATED INTENSITIES CALCULATED BY THREE METHODS 


FOR PROFILES A TO F WITH A HOT WALL a 


Method of calculation 

Temperature distribution 


1 

7 


Integrated intensity, 

I, 

W/ (cm 2 )(steradian) 

Error from 
exact I, 
percent 

Integrated intensity, 

I, 

W /(cm 2 ) (steradian) 

Error from 
exact I, 
percent 


Profile A 




Exact 

One wave -number group* 5 ap^ 
One wave -number group,* 5 one 
origin group a 0 

2.63X10 6 

1.76 

2.22 


-33 

-15 

3.30X10 5 

3.90 

3.55 

18 

8 


Profile B 




Exact 

One wave -number group* 5 ap^ 
One wave -number group,* 5 one 
origin group a 0 

3.91X10 6 

2.60 

3.49 

-34 

-11 

1.73X10 5 
2.45 
1. 68 

42 

-3 


Profile C 




Exact 

One wave -number group* 5 a pl 
One wave -number group,* 5 one 
origin group a e 

4. 36X10 6 
2.49 
4. 18 


-43 

-4 

7. 99X10 4 
2.57X10 5 
7. 80X1 0 4 

222 

-2 


Profile D 




Exact 

One wave -number group* 5 a p ^ 
One wave -number group,* 5 one 
origin group a 0 

1.47X10 6 

1.48 

1.83 

1 

24 

3.22X10 5 

3.48 

3.04 

8 

-6 


Profile E 




Exact 

One wave -number group* 5 a pl 
One wave -number group,* 5 one 
origin group a 0 

1.53X10 6 

1.66 

1.84 

9 

20 

3. 14X10 5 
3. 14 
3.00 

0 

-4 


Profile F 




Exact 

One wave -number group* 5 a pl 
One wave-number group,* 5 one 
origin group a fi 

2.72X10 6 

1.70 

2.03 

-37 

-25 

3.65X10 5 

3.88 

3.79 

6 

4 


a Hot, black wall is at distance L along ray of 50 cm and is at same temperature as adjacent gas. 
^Wave-number group extended 0 to 400 000 cm”*. 
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TABLE IV. - INTEGRATED INTENSITIES CALCULATED BY FIVE 


METHODS FOR PROFILES G AND H WITH A COLD WALL* 
[Temperature distribution designated 8 .] 


Method of calculation 


Exact 


b 


One wave -number group a pi 
Two wave number groups 0 
One wave -number group, b 
origin group a 


origin group a 


Profile H 


Exact 

One wave -number group* 3 a pl 

Three wave -number groups* 1 a p ^ 
One wave-number group, b one 
origin group a 

® d 

Three wave -number groups, 
one origin group a 


5. 09x10 
5.00X10 5 

1.24X10 5 

6.88X10 4 


6. 49x10 



Integrated intensity, 

Error from 



exact I, 


9 

W/(cnO(steradian) 

percent 

Profile G 



4. 17x10® 

— 

i pi 

7.53X10 1 

-100 

a pi 

4. 17x10° 

0 

one 

1.63x10* 

-96 

one 

4. 17X10 5 

0 


881 

143 

35 

27 


a Cold wall is at distance L along ray of 50 cm . 
b Wave-number group extended 0 to 400 000 cm” 1 . 

°Wave-number groups extended 0 to 20 000 and 20 000 to 
400 000 cm' 1 . 

^Wave-number groups extended 0 to 45 830, 45 830 to 75 830, and 
75 830 to 400 000 cm” 1 . 



TABLE V. - INTEGRATED INTENSITIES CALCULATED BY FIVE METHODS FOR BINARY 
MIXTURE A, BINARY MIXTURE B, AND TERNARY MIXTURE WITH A COLD WALL a 


Method of calculation 


Temperature distribution 



8 


9 


Integrated intensity, 

Error from 

Integrated intensity, 

Error from 


, i. 

exact I, 

L 

exact I, 


W/(cm <2 )(steradian) 

percent 

W/(cm^)(steradian) 

percent 


Binary mixture A 



Exact 

2.71X10 5 

— 



One wave -number groups a pl 

1.93 

-29 



Three wave -number groups 0 a pl 

3.32 

23 



Three wave -number groups, c one 

1.89 

-30 



origin group a_ 

e c 

Three wave -number groups, two 

2.59 

-5 



origin group a 0 






Binary mixture B 



Exact 

1. 55X10 5 

— 



One wave -number groups a pl 

7.75 

401 



Three wave -number groups 0 a pl 

3.03 

96 



Three wave -number groups, 0 one 

1.84 

19 



origin group a 

Three wave -number groups, two 

1.98 

28 



origin group a 0 






Ternary mixture 



Exact 



5. 68X10 5 

— 

One wave -number group* 3 a pl 



1.77 

-69 

Three wave -number groups 0 a pi 



6.48 

14 

Three wave -number groups, 0 one 



2.88 

-49 

origin group a_ 

e n 

Three wave-number groups, two 



4.69 

-17 

origin group a g 






a Cold wall is at distance L along ray of 50 cm for binary mixtures and 150 cm for ternary mixture. 
^Wave-number group extended 0 to 400 000 cm'*. 

c Wave-number groups extended 0 to 100 000, 100 000 to 200 000, and 200 000 to 400 000 cm”*. 
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TABLE VI. - LINE POSITIONS AND LINE STRENGTH 


PARAMETERS FOR GASES a AND b OF MIXTURES 


Gas a 

Gas b 

Line position, 

Line strength 

Line position, 

Line strength 

*o> 

parameter, 

* 0 ’ 

parameter, 

cm ” * 

J 

cm”* 

J 

9 994 

1.2022 

33 606 

0. 7986 

19 362 

1. 9862 

56 418 

.2602 

37 438 

2.3076 

73 998 

2.8975 

72 834 

.3478 

82 341 

.4347 

90 324 

.0792 

125 624 

.4061 

119 202 

1.5181 

146 670 

3. 1317 

123 356 

.7789 

159 466 

.4395 

129 578 

.8320 

176 833 

1.3123 

131 880 

. 1995 

177 817 

.8202 

137 408 

.1505 

206 413 

.5250 

150 498 

.3421 

214 598 

1.4836 

182 760 

.7748 

260 527 

.5263 

206 455 

.6532 

274 945 

.7483 

220 791 

. 4057 

282 562 

1.2592 

308 420 

. 1390 

285 207 

.4235 

338 452 

.2008 

291 665 

2.2219 

343 178 

.0245 

309 732 

.0092 

347 055 

.0956 

339 333 

1.4360 

355 507 

.9856 

388 435 

.5223 

366 994 

.0218 

388 579 

2.0375 



TABLE Vn. - LINE POSITIONS AND LINE STRENGTH 


PARAMETERS FOR GASES d AND e OF MIXTURES 


Gas d 

Gas e 

Line position, 

Line strength 

Line position, 

Line strength 

V 

parameter. 

V 

parameter, 

cm"^ 

J 

cm"^ 

J 

5 499 

0. 1767 

69 134 

0. 4401 

44 935 

.5048 

76 387 

.0758 

83 586 

.2458 

104 551 

3.0978 

84 226 

1.7403 

108 506 

.2534 

91 882 

1. 4745 

128 839 

. 1772 

98 932 

.4911 

164 442 

2.9064 

100 973 

3.2213 

195 023 

.7689 

116 644 

1. 1733 

195 235 

.6671 

118 050 

1.4220 

217 783 

1.9388 

125 507 

1. 1864 

218 185 

2. 1993 

128 079 

1. 1195 

218 264 

1.8557 

154 744 

3.2730 

246 373 

.3219 

186 333 

.1222 

272 675 

.9993 

195 654 

2. 1874 

283 060 

.8711 

235 237 

.7292 

306 927 

.4269 

310 601 

.3628 

337 134 

.6405 

321 790 

1.6439 

374 029 

1.5674 

331 851 

1. 1472 

380 010 

.5573 

359 631 

1. 1099 

385 080 

1. 3798 

375 420 

.9530 

398 209 

.9496 
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Figure 5. - Lorentz-line spectral absorption coefficients designated profiles B 
and C. Normalized Planck function is given for reference and is for 25 000° K. 


Normalized Planck function, Brj/B- 





Normalized Planck function, Bc-r/B- 



Wave number, v, cm" 1 

Figure 7. - Spectral absorption coefficient designated profile G. Normalized 
Planck functions for two temperatures are given for reference. 



Figure 8. - Variable-width Lorentz-line spectral absorption coefficient designated 
profile H. Normalized Planck function is given for reference. 
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Ratio of approximate to exact integrated 
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